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It was recently reported that there may exist monochromatic 7-ray emission at ~ 130 GeV from 
the Galactic center in the Fermi Large Area Telescope data, which might be related with dark matter 
(DM) annihilation. In this work we carry out a comprehensive check of consistency of the results 
with the DM annihilation scenario, using the Fermi observations of the inner Galaxy, Galactic halo, 
clusters of galaxies and dwarf galaxies. The possible effects on DM annihilation of substructures 
based on recent simulations are also taken into account. We find that: 1) the constraints on the 
annihilation cross section of continuous 7-ray emission from the Galactic center are as stringent as 
the "natural" scale assuming thermal freeze-out of DM; 2) the Fermi data of Galactic halo seems 
to have tension with the DM annihilation interpretation of the ~ 130 GeV emission from the inner 
Galaxy; 3) possible concentration of photons in 120 — 140 GeV from nearby clusters is revealed; 
4) constraints from galaxy cluster (with substructures) are marginally consistent with the DM 
annihilation scenario to explain the ~ 130 GeV emission, and the constraints from dwarf galaxies 
are weaker. 

PACS numbers: 95.35.+d,95.85.Pw 



I. INTRODUCTION 

Through analyzing the public Fermi Large Area Tele- 
scope (Fermi-LAT) data, several groups reported the 
hints of monochromatic 7-ray emission with energy E 
130 GeV [H-1H. The morphology of the potential line 
emission is still in debate. In some works the results 
showed that there were various regions lying basically in 
the Galactic plane having excess in 120 — 140 GeV 0,3, 
while in another work the excess concentrated in the in- 
ner 5° of the Milky Way center [j| . The results stimulated 
active discussions of the possible dark matter (DM) ori- 
gin . Alternatively the astrophysical explanations of 
the line- like 7-ray emission were also proposed 0, [l(| ■ 

Due to the importance of the high energy line emission, 
many works were trying to test the line emission (and 
its DM interpretation) with Fermi-LAT observations of 
continuous 7-rays [HHHI, the line emission from dwarf 
galaxies and clusters of galaxies [HI, the unassociated 
Fermi sources [l7j . and the future detection with high 
energy resolution telescope The basic results show 
that the allowed cross section of the DM annihilation final 
states giving rise to the continuous emission can only be 
about O(10 2 ) times larger than that of the line emission 
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, which needs to be considered when constructing 
the DM models. The recent search for line emission in 
the Milky Way halo by Fermi collaboration with two- 
year PASS 6 data showed no indication of signals fl9| . 
and the upper limits seemed to be marginally in conflict 
with the results found in 0, H| • Furthermore, no signal 
from dwarf galaxies was found and the constraints on DM 
annihilation into monochromatic 7-rays were consistent 
with the results found in the inner Galaxy [15| . For the 



galaxy clusters, however, it was claimed a ~ 3er signal of 
possible 130 GeV line emission by Hektor ct al. [la ]. 

In this work we try to make a comprehensive test of 
the DM scenario of the line emission, through analyzing 
the Fermi-LAT 7-ray data in the Galactic center, the 
Milky Way halo, dwarf galaxies and galaxy clusters. In 
our analysis the information of the spatial distribution of 
DM is taken into account. We further include the effects 
of substructures on the DM annihilation in the analysis, 
both the enhancement of annihilation luminosity and the 
change of surface brightness distribution pfjj ]. 

This paper is organized as follows. In Sec. II we will 
discuss the constraints on the continuous 7-ray emission 
from the Galactic center region. In Sec. Ill, IV and 
V we give the results of the line search from the Milky 
Way halo, galaxy clusters and dwarf galaxies respectively. 
Finally Sec. VI is our conclusion and discussion. 



II. CONTINUOUS GAMMA-RAY EMISSION 
FROM THE INNER GALAXY 

As DM is a neutral particle it can not couple to photons 
directly. Therefore DM annihilate into to monochromatic 
photons through loop processes with charged virtual par- 
ticles in the loop. If DM annihilate into the charged par- 
ticles at the tree level, such as quarks, charged leptons 
or gauge bosons, these annihilation products may induce 
significant continuous 7-ray flux. Using the Fermi data 
we can set a strong constraint on the DM annihilation 
cross section into these final states. 

The current most stringent limits on the DM annihi- 
lation cross section set by the Fermi collaboration were 
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derived from the analysis of ten dwarf satellite galaxies 
with 24 months of Fermi data [2lJ . The Galactic center 
(GC) has been thought to be the best target for probing 
the DM annihilation signals. It may have larger DM an- 
nihilation rate than dwarf galaxies, and is a good target 
for searching for the continuous 7-rays induced by DM. 
In this section, we will use the Fermi data from the in- 
ner Galaxy to place limits on the DM annihilation cross 
section. This analysis is important for building the DM 
model to explain the monochromatic 7-ray emission as it 
has to satisfy the continuous 7-ray limit simultaneously. 

In this analysis we use the 3.7 years Fermi-LAT data 1 
recorded from 4 August 2008 to 18 April 2012, with 
the Pass 7 photon selection. The "SOURCE" (ev- 
class=2) event class is selected, and the recommended 
filter cut " (DATA_QUAL== 1 ) && (LAT_CONFIG==l) 
&& ABS(ROCK_ANGLE)< 52" is applied. The en- 
ergy range of events is restricted from 500 MeV to 300 
GcV, and the region-of-interest (ROI) is adopted to be 
a 10° x 10° box centered around the GC. We carry out 
the binned likelihood analysis with the LAT Scientific 
Tools v9r23pl. The instrument response function used is 
"P7S0URCE_V6" . For the diffuse background, we use the 
Galactic diffuse model gal_2yearp7v6_v0 . f its and the 
isotropic diffuse spectrum iso_p7v6source .txt provided 
by the Fermi Science Support Center 2 . 

The data are binned with 30 energies bins logarithmi- 
cally spaced and 100 x 100 spatial bins with size 0.1° each. 
For the likelihood analysis we include the point sources 
within 5° around the Galactic center based on the 2-year 
LAT source catalog (using the user-contributed software 
make2FGLxml.py) [22(. For other sources located in the 
ROI we fix their parameters to be the values of the LAT 
source catalog. Before taking the DM contribution into 
account, we first make a global fit to derive the spectral 
parameters of the point sources. All the spectral param- 
eters of the sources within 5° of the Galactic center and 
the normalizations of the diffuse backgrounds are left free 
during the fit. Then we add the DM component as a dif- 
fuse source, and re-do the fit to get the upper limits of the 
DM contribution. The free parameters of the latter fit 
include the normalizations of the point sources within 5° 
of the Galactic center, the normalizations of the diffuse 
backgrounds and the DM component. 

For the DM component we discuss W + W~ , bb, p + p~ 
and t + t~ final states. The photon production spectra 
are calculated with the P YTHIA simulation package [23| . 
The spatial density distribution of DM is assumed to be 
either an Navarro- Frenk- White (NFW, HI) profile 
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with r s = 20 kpc and p s — 0.35 GeV cm 3 , or an Einasto 
profile H| 



p(r) = p s exp 
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where a = 0.17, r s = 20 kpc and p s = 0.08 GeV cm~ 3 . 
The scale density of both theprofiles is to give a local 
DM density ~ 0.4 GeV cm" 3 (2f3|. 

In the analysis we also take the DM substructures into 
account. We employ the results from Aquarius simula- 
tions [23] and the fitted annihilation rate from the sub- 
structures in Ref. 12811 
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where M max « 10 M is the maximum mass of sub- 
haloes in the Milky Way halo, M Tes ps 3 x 10 5 M is the 
resolution mass of the simulation, and the minimum sub- 
halo mass is assumed to be M m j n ps 10~ 6 Mq. For the 
inner Galaxy discussed in this section, the substructures 
play little role in the enhancement of the annihilation 
signal. 

We project the density square into a 2-dimcnsional sur- 
face map to give the spatial template of the DM annihila- 
tion. With the spatial template and the energy spectrum, 
we use the above described likelihood analysis method to 
derive the upper limits of the DM component. The only 
parameter of the DM component is the normalization 
factor. 
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FIG. 1: Constraints on the DM annihilation cross sections to 
W + W~ (solid), bb (dashed), p + p~ (short-dashed) and t + t~ 
(dash-dotted), for 10° x 10° region around the Galactic cen- 
ter. Lines are for NFW profile, and lines with points are for 
Einasto profile. 



1 http:/ /fermi. gsfc.nasa.gov/ssc/data 

2 http:/ /fermi. gsfc.nasa.gov/ssc/data/access/lat/Background- 
Modcls.html 



The results of the constraints are shown in Fig. Q] It 
can be seen that for DM particle with mass less than 
hundreds of GeV most of the constraints can reach the 
"natural" scale, 3 x 10~ 26 cm 3 s _1 , assuming DM is the 
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thermal relic of the early universe. These constraints are 
so far the most stringent constraints on the DM annihi- 
lation with Fermi-LAT 7-ray data. Compared with fl9| . 
our constraints are more stringent because we use the 
likelihood fit involving the astrophysical background, in- 
stead of assuming the DM component not to exceed the 
observational data conservatively. Also in [29| the con- 
straints were given directly by comparing the DM con- 
tribution with the observed excess in the Galactic center 
region with the known sources and diffuse emission sub- 
tracted. 

For DM mass - 130 GeV and NFW (Einasto) pro- 
file, the constraints on the cross sections of W + W~ , 
bb, n+yT and t+t" arc 7.4 x 10" 26 (3.7 x 1CT 26 ), 

5.1 x 10~ 26 (2.6 x 1CT 26 ), 4.8 x 10" 26 (2.0 x 10" 26 ) and 
3.6 x 10~ 27 (1.6 x 10~ 27 ) cm 3 s" 1 respectively. Com- 
pared with the cross section to monochromatic 7-rays 
ov mono « 2.3 x 10~ 27 (1.3 x 10~ 27 ) cm 3 s" 1 0, the con- 
tinuous cross section can only be larger by a factor of 

1.2 — 30. Our results are slightly stronger than that de- 
rived in fl2T4l4l | , which may be due to the use of spatial 
information of DM distribution in this analysis. 

Such strong constraints on the DM annihilation cross 
section of the tree level processes are very impor- 
tant when constructing the DM models to explain the 
monochromatic emission. In the following we discuss the 
implications of the constraints briefly 

As we have mentioned, DM annihilate into monochro- 
matic photons through higher order processes with 
charged particles running in the loop. These charged 
particles will induce large continues 7-ray flux. Here we 
define a critcrium as the ratio of the tree level process 
into continuous 7-rays and the loop process into line 7- 
ray, R = <TV CO nt/vv mono 

= av xxl av ii where X denote 
a charged particle. To satisfy the Fermi constraints, R 
should be smaller than 1.2 — 30, depending on different 
annihilation modes. 

The charged particle, X, in the loop can be a SM par- 
ticle and/or an exotic particle in new physics. Since 
the loop processes are generally suppressed, the large 
DM annihilation cross section into monochromatic pho- 
tons ~ 10~ 27 cm 3 s _1 often requires additional enhance- 
ment mechanism. If all the charged particles in the loop 
are SM particles, an additional narrow resonance medi- 
ator A is useful to enhance DM annihilation cross sec- 
tion [[|, [llj . Using the narrow width approximation, R 
can be determined by the branching ratios of A decay 
BR(A -» XX)/BR{A 77). In general, R is pro- 
portional to 1/a 2 due to the loop factor and two QED 
couplings. If X are quarks or leptons, R would be en- 
hanced by an additional factor of m^/m^. In this case, 
the continuous 7-ray flux is much larger than monochro- 
matic flux, for instance by a factor of > 10 6 , and are 
certainly excluded by Fermi results. If X is W boson, 
R is ~ O(10 4 ) and can be excluded by continuous 7-ray 
observations too. 

If the mediated charged particles in the loop include 
both SM particles and exotic particles in new physics, 



the estimation of R depends on the details of model, 
and is more complicated. A familiar example is super- 
symmetry in which the neutralinos can annihilate into 
photons through box loop including chargino-IU boson 
and sfermion-fermion contributions (30|. If the main 
component of DM is Wino or Higgsino, the diagram in- 
cluding chargino-IU boson contributions would be dom- 
inated. In such cases, the large Wino/Higgino-chargino- 
W couplings are also helpful to obtain required anni- 
hilation cross section. However, R is still as large as 
O(10 2 ) - O(10 3 ) at least [3, and can be excluded by 
Fermi continuous 7-ray limits. 

Therefore, these constraints suggest the charged parti- 
cles in the loop have to be heavier than the DM particle 
so that the tree level process into the charged particles is 
forbidden. For example, a chargino triangle loop can be 
used to explain required DM annihilation cross section 
into monochromatic photons 0, H| . In this case, the pro- 
cesses producing line and continuum 7-ray spectra are 
independent. Therefore, no excess continuous 7-ray flux 
would be induced by adjusting model parameters. 

Finally, we give another comment on the importance 
of the continuous 7-ray constraints. Any successful DM 
model needs to explain the observed DM relic density 
Q,h 2 ~ 0.11. If DM is thermally produced through 
velocity independent annihilation process in the early 
universe, the cross section is about the "natural" scale 
3 x 10~ 26 cm 3 s _1 . If the future upper-limits derived 
from continuous 7-ray observations on DM annihilation 
cross section are lower than the "natural" scale, it means 
the DM should be totally or partially produced through 
other mechanisms. For example, the corrected DM relic 
density can be obtained by co-annihilation, or velocity 
suppressed annihilation jll| . or non-thermal production 
mechanism (2^|. These results will be important for the 
DM model building and DM detections. 



III. LINE EMISSION FROM THE GALACTIC 
HALO 

The tentative 7-ray line emission was found mainly in 
the inner Galaxy, which is expected in the DM annihi- 
lation scenario. As a consistency check, the search for 
the line emission in the Galactic halo should be of great 
importance, especially in case that there might be sig- 
nificant contribution from DM substructures as shown 
by the cold DM simulations. We use the same period 
and data version as in Sec. II, but the whole sky re- 
gion is adopted. Both the event classes "SOURCE" and 
"CLEAN" are analyzed. Since we focus on the high 
energy line emission, the energy range in this analysis 
is adopted to be [20, 300] GeV. We choose two sky re- 
gions in the analysis. The first one (hereafter Reg. 1) is 
|6| > 10°, which can exclude most of the Galactic plane 
and keep most of the halo region. The second one (Reg. 
2) is |6| > 10° plus \l\ < 10°, \b\ < 10°, due to the fact 
that DM will concentrate in the central Galaxy and such 
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an adoption of the sky region may keep more potential 
DM signal. The latter adoption is the same as that in 
Ref. . We also test the effect of point sources, through 
masking 1° radius of each point source based on the sec- 
ond Fermi-LAT source catalog (22[. 

To search for the 7-ray line emission, we perform a 
likelihood fit to the fluxes using a power-law background 
together with a monochromatic 7-ray line. The energy 
resolution of Fermi-LAT detector (as given in [l9[ with 
an average of the angular acceptance) has been taken 
into account in the fit. Through the fit we can derive the 
95% confidence level upper limits on the line emission, 
which can be translated into constraints on the DM an- 
nihilation cross section to monochromatic photons. To 
get the constraints on the cross section, we first compute 
the dimcnsionlcss astrophysical J-factor, defined as 
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where p© = 0.4 GeV cm~ 3 is the local density, Rq = 8.5 
kpc is the distance from the Galactic center to the Earth, 
and l.o. s. means integral along the line of sight. The J- 
factors for NFW and Einasto profiles, with and without 
substructures are listed in Tabic HI The numbers in paren- 
thesis are the results after masking 1° radius around each 
point source of the 2-ycar LAT source catalog p2j . 

Constraints on the DM annihilation cross section are 
shown in Fig. O The left panel is for Reg. 1 and the 
right panel is for Reg. 2. Here we adopt the Einasto 
density profile, and the possible effects from substruc- 
tures are not considered. The cross section derived in 
Weniger (2012) is also shown for comparison. The re- 
sults show that it is consistent with Wcniger's result when 
we choose "SOURCE" data as he did in [3, but for the 
"CLEAN" data result there might be tension between 
the inner Galaxy and the Milky Way halo observations 
of the line-like emission if explaining it with DM anni- 
hilation. For NFW profile, the tension is even stronger. 
If the DM substructures are taken into account, the con- 
straints from Milky Way halo will be stronger by a fac- 
tor of ~ 2. However, for the "SOURCE" data, even the 
substructures are included the results are still concictent 
with that by Weniger (2012). The "CLEAN" results of 
Reg. 2 are also consistent with that derived by Fermi- 
LAT collaboration [l9j], in which two-year Pass 6 "UL- 
TRACLEAN" data are used. 



IV. LINE EMISSION FROM GALAXY 
CLUSTERS 

Galaxy clusters (this section) and dwarf galaxies (next 
section) are also potential good targets for probing DM 
signals. In this section we analyze the Fermi-LAT data 
on a series of nearby clusters to search for the poten- 
tial line emission. The data version and time window 
are the same as in the above discussions, but the en- 
ergy range is set from 200 McV to 300 GeV. The ROI is 



adopted to be 14° x 14° box centered in the target galaxy 
cluster. Seven clusters, Fornax, AWM7, M49, NGC4636, 
Centaurus, Coma and Virgo are choosen as our targets. 
The mass of Virgo is adopted from [3l[ , and for others 
the masses are adopted from the extended HIFLUGCS 
X-ray catalog [H, HH . 

Fig. |3] shows the individual photons with arrival direc- 
tions within the virial radius of these clusters. The left 
panel is for "SOURCE" events, and the right panel is for 
"CLEAN" events. It is interesting to find that for the 
"CLEAN" events, photons with energies between 120 and 
140 GeV are more abundant than nearby energy ranges. 
We checked for the "ULTRACLEAN" events and find 
identical results with the "CLEAN" events. From Fig. 
|3] it seems hard to claim the existence of line emission. 
Recently it was found a 3er evidence for the existence of 
~ 130 GeV line emission from galaxy clusters in fl6j |. 

We then derive the constraints on the DM annihilation 
cross section to 7-ray line from the galaxy cluster obser- 
vations. The DM signals from galaxy clusters are mod- 
eled as extended sources as in 34 1. with surface bright- 
ness 



S sm (9) ex 



p 2 am (l,9)dl, 
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where 9 is the angle between l.o.s. and the cluster cen- 
ter, p sm is the density profile of the smooth halo, which 
is assumed to be NFW profile. The mass-concentration 
relation is adopted to be (35| 
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The DM substructures may play a significant role on 
the enhancement of annihilation signal from clusters [20L 
l36| | . Following Ref. [2(3| we adopt the boost factor of 
substructures as a function of cluster virial mass Af 2 oo 

KMaoo) = 1-6 x 10- 3 (A/ 200 /M Q ) a39 , (7) 

and the surface brightness profile of subhalo emission 

16&(M 2 oo)S'nfw 1 
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where S"nfw is the integral flux of the smooth NFW halo 
of the cluster, cIa is the angular diameter distance. 

To study the 7-ray line, we need to take into account 
the energy resolution of the detector. We use an energy 
dependent resolution from (l9| , where the energy resolu- 
tion was derived through Monte-Carlo simulations with 
an integral over the angular acceptance. Then we use 
the binned likelihood analysis tool as in Sec. II to de- 
rive the upper limits of DM annihilation final state 77 . 
The point sources of the second Fermi-LAT catalog [22( 
located in the ROI are modeled simultaneously, with the 
normalizations free. The normalizations of the diffuse 
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TABLE I: Jaq of the DM annihilation in the Milky Way halo. In parenthesis are the results after masking point sources. 
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NFW+sub 


EIN EIN+sub 
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Reg. 2 


20.2(17.7) 
31.4(22.5) 


43.9(38.7) 
55.4(43.6) 
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42.2(27.8) 66.2(48.9) 
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FIG. 2: Constrains on the DM annihilation cross section to monochromatic 7-rays. Left panel is for Reg. 1 and right panel is 
for Reg. 2. The DM density profile is assumed to be Einasto profile, and no substructures are included. 
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FIG. 3: Scattering plot of photon energies and deviation angles from the cluster center, for the events within virial radius of 
each cluster. Left panel is for "SOURCE" events, and right panel is for "CLEAN" events. 



backgrounds are also left free in the likelihood fit. Be- 
sides the constraint from individual galaxy clusters, we 
also do a combined likelihood analysis of all the seven 
galaxy clusters simultaneously to combine the statistical 
power of these different target regions [34| • Constraints 
on the DM annihilation cross section are shown in Fig. 
31 for the cases without (left) and with (right) substruc- 
tures. It can be seen that the constraints from galaxy 
clusters are quite weak and are still consistent with the 
DM interpretation of the inner Galaxy line emission, even 
for the case with significant substructure boost. 



V. LINE EMISSION FROM DWARF GALAXIES 

Same as the analysis of galaxy clusters, we perform the 
search for line emission with Fermi-LAT data of dwarf 
galaxies in this section. In Ref. [ll| the authors have 
done a similar analysis and found no signal of line emis- 
sion from the dwarf galaxies. Here we redo the analysis 
with a slightly different way 3 from that of Ref. [ll| , for 
the seek of completeness and independent check. We 
use the same data selection conditions as above in the 



Specifically, our method is the same as the analysis of galaxy 
clusters, with point sources in the ROI and diffuse backgrounds 
involved. 
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FIG. 4: Constraints on DM annihilation cross section to 77 from clusters without (left) and with (right) substructure enhance- 
ment. 



galaxy clusters analysis. The targets chosen are Bootes 
I, Draco, Fornax, Sculptor, Segue 1, Sextans and Ursa 
Minor, which are the same as in Ref. (l5j . For each 
dwarf galaxy, we model it as point source and the anni- 
hilation J-factor is adopted to be the best fit value given 
in [33]. Then we perform the binned likelihood analysis 
to search for the line emission and derive the constraints 
on monochromatic 7-ray line flux. 

We have found no significant 7-ray line emission in the 
target regions. The constraints on (av) xx _,. 77 are shown 
in Fig. [5] We can see that the upper limits are still far- 
away from the best fit point from Ref. Q. This result is 
consistent with the upper limit derived in fl5j , though we 
use a different method to deal with background photons. 
At the present time the constraint from dwarf galaxies 
can neither confirm nor exclude the DM-induccd line in- 
terpretation of the Galactic center data. 
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FIG. 5: Constraints on DM annihilation cross section to 77 
from dwarf galaxies. 



VI. CONCLUSION AND DISCUSSION 

The recently reported tentative 7-ray line emission 
from the Fermi-LAT observations in the inner Galaxy p]- 
[f| has invoked many discussions of the DM annihilation 
signals. In this paper we do a comprehensive analysis 
using the Fermi-LAT data in the Galactic center region, 
Galactic halo, galaxy clusters and dwarf galaxies, to test 
the DM annihilation interpretation of this line emission. 

Using the Fermi-LAT data in the Galactic center re- 
gion, we get strong constraints on the continuous 7-ray 
emission of DM annihilation to W + W~, bb, fi + fJ,~ and 
t + t~ final states, which set useful constraints on DM 
models to explain the 130 GeV line emission. 

We further perform the search for line emission with 
Fermi-LAT data in the Milky Way halo. The results seem 
to have tension with the DM annihilation interpretation 
of the 130 GeV line emission if the DM density profile fol- 
lows Einasto or NFW profiles, if the "CLEAN" events are 
employed. Considering the possible enhancement effect 
due to DM substructures in the halo, the tension is even 
more stringent. To consistently understand the results in 
the inner Galaxy and the Galactic halo within DM anni- 
hilation scenario, one may need to assume cuspier den- 
sity profile of DM in the Galactic center (e.g., accreted 
by the central supermassive black hole [231). However, 
the "SOURCE" events are still consistent with the DM 
scenario to explain the 130 GeV line emission in Ref. @]. 

Galaxy clusters and dwarf galaxies are also studied to 
further test the DM interpretation of the 130 GeV line 
emission. It is interesting to find that there is a con- 
centration of photons in 120 — 140 GeV from the six 
nearby clusters (Fig. 4) for the "CLEAN" and "ULTR- 
ACLEAN" classes of events. However, the result seems 
hard to indicate an excess of line emission. Therefore 
we set upper limits on the DM annihilation cross section 
to monochromatic 7-ray line, which are consistent with 
the claimed results in the inner Galaxy 0, even when t he 
large boost factor of substructures are taken into account. 
No signal of line emission from dwarf galaxies is found, 
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and the constraints on DM annihilation cross section to 
monochromatic 7-ray line are derived. Compared with 
the galaxy clusters, the constraints are weaker for dwarf 
galaxies, and are consistent with the claimed results in 
the inner Galaxy @. 

Acknowledgements 

We thank Yi-Zhong Fan and Christoph Wenigcr 
for helpful discussion and useful comments on the 



manuscript. This work is supported by National Natural 
Science Foundation of China under grant Nos. 11075169, 
11073024, 11105155, 11105157, the 973 project under 
grant No. 2010CB833000, and Chinese Academy of Sci- 
ences under grant No. KJCX2-EW-W01. XLC acknowl- 
edges the support from John Templeton Foundation & 
NAOC Beyond the Horizon program. 



[1] T. Bringmann, X. Huang, A. Ibarra, S. Vogl and C. 
Weniger, larXiv:1203.13T2l 

[2] C. Weniger, larXiv: 1204.2797] 

[3] E. Tempel, A. Hektor and M. Raidal, larXiv:1205.T045l 
[4] A. Boyarsky, D. Malyshev and O. Ruchayskiy, 

larXiv:1205.4700l 

M. Su and D. P. Finkbeiner. laxXiv: 12061616 



[5] 
[6] 



['] 

[8 
[9 
[10 

[11 

[12 

[13 
[14 
[15 

[16 

[17 
[18 

[19] 



M. 
O. 

-C. 
and 



Cline, 
Seto, 
Park, 
W.-I. 



E. Dudas et al, larXiv: 1205. 1520 J. 
larXiv:1205.2688l K.-Y. Choi and 
larXTv : 1205.3276 B. Kyae and J, 
arXiv: 1205.415 1] H. M. Lee, M. Park, 
Park, arXiv:120 O675l A. Rajaram an, T. M. P. Tait, 
and D. Whiteson, larXiv:1205.4723l B. Samir Acharya 
et aL larXiv:1205.5789l X. Chu et al, arXi v:1206.2"279l 
L. Feng, Q. Yuan and Y. Z. Fan, [a rXiv: 120 6^2910] 
J. H. Heo and C. S. Kim, larXiv: 1207. 1341] I. Oda, 
larXiv:1207.1537l R. Z. Y ang et al, larXiv:1207.1621l 
Mads T. Frandsen et al, larXiv: 1 207.3971 ; Jong-Chul 
Park and Seong Chan Park, larXiv:1207.498T1 

D. Das, U. Ellwanger and P. Mitropoulos, 

larXiv:1206.2639l 

Z. Kang et qZ., larXiv:1206.2863l 

S. Profumo and T. Linden. larXiv:1204.6047l 

F. Aharonian, D. Khangulyan, and D. Malyshev, 
larXiv:1207.0458l 

M. R. Buckley and D. Hooper, larXiv: 1205.681T1 
W. Buchmuller and M. Garny, larXiv: 1 206.7056 

T. Cohen et al, larXiv:1207.0800l 

I. Cholis, M. Tavakoli, P. Ullio, larXiv:1207.1468l 

A. Geringer-Sameth and S. M. Koushiappas, 

larXiv:1206.0796l 

A. Hektor, M. Raidal and E. Tempel, larXiv:1207.446"6l 

M. Su and D. P. Fi nkbeiner, larXiv:1207. 7060] 

Y. Li and Q. Yua n, larXrv: 1206.2241] L. Bergstrom et al, 

larXiv:1207.6773l 

Fermi LAT Collaboration, M. Ackermann et. al, 
larXiv:1205.2739l 



[20] L. Gao et. al, Mon. Not. Roy. Astro n. Soc. 419 (2012) 
1721; L. Gao et. al. larXiv:1201. 19401 

[21] Fermi LAT Collaboration, M. Ackermann et. al, Phys. 
Rev. Lett. 107 (2011) 241302 

[22] Fermi LAT Collaboration, P. L. Nolan et. al, Astrophys. 
J. Supp. 199 (2012) 31 

[23] T. Sjostrand, S. Mrenna and P. Skands, J. High Energy 
Phys. 5 (2006) 26 

[24] J. F. Navarro, C. S. Frenk and S. D. M. White, Astro- 
phys. J. 490 (1997) 493 

[25] J. Einasto, Trudy Astrofizicheskogo Instituta Alma-Ata 
5 (1965) 87 

[26] R. Catena and P. Ullio, J. Cosmol. Astropart. Phys. 08 
(2010) 004 

[27] V. Springel et. al, Mon. Not. Roy. Astron. Soc. 391 

(2008) 1685 

[28] Q. Yuan et. al. larXiv:1203.5636l 

[29] D. Hooper and T. Linden, Phys. Rev. D 84 (2011) 123005 

[30] L. Bergstrom and J. Kaplan, Astropart. Phys. 2 (1994) 
261; G. Jungman and M. kamionkowski, Phys. Rev. D 
51 (1995) 3121; Z. Bern, P. Gondolo and M. Perelstein, 
Phys. lett. B. 411 (1997) 86; L. Bergstrom and P. Ullio, 
Nucl. Phys. B. 504 (1997) 27 

[31] A. Pinzke, C. Pfrommer and L. Bergstrom, Phys. Rev. 
D 84 (2011) 123509 

[32] T. H. Reiprich and H. Boehringer, Astrophys. J. 567 
(2002) 716 

[33] Y. Chen et. al, Astron. Astrophys. 466 (2007) 805 

[34] X. Huang, G. Vertongen and C. Weniger, J. Cosmol. As- 
tropart. Phys. 01 (2012) 042 

[35] A. R. Duffy et. al, Mon. Not. Roy. Astron. Soc. 390 
(2008) L64 

[36] J. Han et. al. larXiv:1201.1003l 

[37] Fermi LAT Collaboration, M. Ackermann et. al, Phys. 

Rev. Lett. 107 2011 241302 
[38] P. Gondolo and J. Silk, Phys. Rev. Lett. 83 (1999) 1719 



